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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

INVESTIGATION OF THE LIFT, CENTER OF PRESSURE, AND 
DRAG OF A PROJECTILE AT A MACH NUMBER OF 8.6 
AND A REYNOLDS NUMBER OF IT MILLION 
By Thomas N. Canning 

SUMMARY 


Models of a high-speed projectile configuration were tested in free 
flight to determine the lift-curve slope, center of pressvire, and drag 
at Mach numbers near 8.6 and at a Reynolds number of 17 million. These 
results were compared with predictions based on available theory and 
experimental results. The drag of a similar cone model was also measured. 


INTRODUCTION 


The aerodynamic consequences of improving the static longitudinal 
stability of a 5° half -angle cone by adding a short spike of heavy mate- 
rial onto the nose are considered in this paper. The resulting shape 
(see fig. l) was enough different from a cone that some changes in aero- 
dynamic properties were expected. The lift on the added cylinder was 
expected to increase the lift at angle of attack and to cancel some of 
the gain in stability by moving the center of pressure forward. It was 
also expected that there would be at least a small region of quasi two- 
dimensional flow just behind the spike which might alter the lift, center 
of pressure, and drag. The presence of the shock wave at the spike base 
was expected to promote boundary- layer transition. The net effects of 
these changes were of sufficient interest that the present experimental 
investigation was undertaken. Tests were made to determine lift, drag, 
and center of pressure of the projectile. The drag of immodified cones 
was also measured for comparison. 

Since the speeds of interest for this body were in the high super- 
sonic Mach number range, the tests were run at the highest Mach number 
compatible with model strength considerations. The tests were made at 
speeds near M = 8.6 in the Ames supersonic free-f light wind tunnel at a 
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Reynolds member of over 17 million, based on model length. In this facil 
ity the models are fired upstream through the test section of an M = 2,0 
supersonic wind tunnel to attain high relative airspeed. 


SYMBOLS 


a 


CD 

^Do 

% 


projection in the x,z plane of the instantaneous accelera- 
tion of model centel* of gravity normal to the tunnel center 
line, ft/sec^ 

drag coefficient, 

oqo 

drag coefficient at a = 0 

lift coefficient corresponding to mean vector sum of a and 3 
in test interval 

dCL 

lift-curve slope, per radian 


Cm 


D 


pitching-moment coefficient about model center of gravity, 
pitching moment 

Slqo 


pitching-moment-ciirve slope 


dCm 

da 


, per radian 


drag, lb 


E,F constants defining variation of a with time 

e base of Naperian logarithms 

f frequency of pitching motion, cps 

I principal moment of inertia about lateral axis through model 

center of gravity, ft-lb sec^ 

k damping constant, sec ^ 

L lift, lb 

I model length, ft 

M test Mach number 


m model mass, slugs 
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Q double integral with respect to time of the angle- of- attack 

function, sec^ 

free-stream dynamic pressure, Ib/sq ft 

Reynolds number based on model length and free-stream conditions 
frontal area of model, sq ft 
time, sec 

projection in x-z plane of the instanteineous velocity of model 
center of gravity normal to tunnel center line, ft /sec 

^g#^cp distance from model nose to center of gravity and center of 
pressure, respectively, ft 

x,y,z orthogonal coordinate system using tunnel center line as x 

axis, ft 

a angle of attack of model relative to local flight path, radians 

Or mean value of vector sum of a and p during test interval, 

radians 

P angle of sideslip of model relative to local flight path, 

radians 

^ 2itf, radl6ins/sec 


Subscripts 


1,2, 3, 4 stations in wind tunnel 


APPARATUS AND TECHNIQUES 
Wind Tunnel 


The tests were conducted in the Ames supersonic free-flight wind 
txmnel. This facility is a short ballistic range inside a variable- 
pressure, supersonic, blowdown wind tunnel. While the tunnel was operat- 
ing at a Mach number of 2.0, the models were fired upstream through the 
15 -foot-long test section at 5600 feet per second from a gun located in 
the diffuser. The aerodynamic data were obtained from a history of the 
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model motion, as recorded by seven shadowgraph stations (four parallel 
to the x,z plane and three parallel to the x,y plane) and a chrono- 
graph. The data from the fo\ir-station group are used to study the motion 
of the models as if the model motion were confined to the x,z plane. 
Details of the data analysis and a discussion of measures taJcen to account 
for motion outside the x, z plane are included in the section entitled 
"Data Reduction." Details of tionnel operation are given in reference 1. 


Models 


The models (fig. l) were machined from 75 ST-6 aluminum bar stock. 
For the projectile the basic 5° half -angle cone was altered by adding a 
cylinder of one third the length, and of 0.19 the base diameter of the 
original cone. The cylinder and nose cone would pres\jmably be made of 
heavy material in a full-size missile to promote stability. The models 
for the present tests were bored out at the base to provide the same 
effect. A conical hole with roxinded bottom was bored in each projectile 
model to place the center of gravity at about the 71-percent point aft of 
the nose. This resulted in a wall thickness of about 0.035 inch, which 
permitted the model to be accelerated about 400,000 gravities in the gun 
without deformation. The nose cones were finished carefully and the 
blmtest one fired was 0.003 inch flat at the tip; most tips were less 
than 0.001 inch flat. 


The surface finish was produced by a fine finishing cut and a little 
hand sanding in the lathe. No effort was made to improve the surfaces 
further and they were satin-like in appearance. Several model surfaces 
were examined in order to evaluate the equivalent roughness at large 
scale. A photograph of a typical model is given in figure 2, and photo- 
micrographs of a model surface and profile are given in figure 3. The 
pictures are representative of the worst conditions noted; for exanple, 
the dip, which is indicated in the profile picture is 0.0001 inch deep by 
0.0030 inch long. This profile is believed to be representative of the 
test models, but no extensive examination of model surfaces was made. 

The models were fired from a high performance .59 caliber gun and were 
held on the sabot by evacuating the space in the base of the model as 
indicated in figure 1. The high-pressure powder gases in the gun chamber 
penetrated the plastic-tape seal, and a small amount of gas flowed into 
the evacuated model. The resistance to acceleration held the model firmly 
in place until the sabot was free of the gun. At this moment the powder- 
gas pressure in the model forced separation from the sabot face. 
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Data Reduction 


The lift-curve slope, center of pressiire, and drag were calculated 
from the record of model motion as a function of time. The component of 
the resultant aerodynamic force parallel to the z axis, that is, the 
lift, was obtained by measuring the acceleration along this direction; 
the center of pressure was foitnd from the natural pitching frequency and 
the lift-curve slope; and the drag was determined by measuring the axial 
deceleration. Details of the lift and center-of -pressure calculations are 
presented below; the details of drag calculation, which is more readily 
visualized, are given in reference 1. For the pvirpose of data reduction, 
linear aerodynamic characteristics are assumed. The effect of nonline- 
arity is treated approximately. 

As mentioned previously -the data from the foxar shadowgraph stations 
recording the model motion in the x,z plane were used as if the model 
were free to move in this plane only. The instantaneous acceleration of 
the model projected, on the x,z plane may be written 


d^z ^ ^LgOgoS 
dt^ m 


which is Newton's law of motion in the x,z plane. This assumes a linear 
lift curve and neglects the contribution of the small lift force due to 
pitching and plunging. 


The time variation of angle of attack in equation (l) can be obtained 
using the angle-of -attack measurements from the shadowgraph pictures and 
two assumptions which define the form of the motion. The assumptions are 
that the restoring moment and damping moments are proportional to the angle 
of attack and pitching rate, respectively. This leads to a variation of 
angle of attack given by 


a 


-kt 

e 


(E cos wt - F sin wt) 


(2) 


Equation (2) is fitted to the observed variation of a with respect to 
time by a least-squares procedure described in reference 2. In this way 
the four unknowns, u, k, E, and F, are evaluated. 


Combining equations (l) and (2) gives the plunging acceleration as 
a function of time 


d^z _ ^ 
dt^ ” dt 


CLcc^IoS 

e 

m 


(E cos wt .- F sin uit) 


( 3 ) 
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Integrating vitli respect to t gives the following equation for vertical 
velocity, w 


e (E cos wt - F sin cot) dt (4) 

in which Wi is the component of the model velocity along the z axis 
at the first station. A second integration gives the equation of verti- 
cal position relative to the initial point zi. The integrals are evalu- 
ated between limits corresponding to the times and positions in two shad- 
owgraph stations . 


dz 


Ci^qoS 


m 


Z 2 = zi + wi (ta - ti) + P ^ (E cos u)t - F sin tut) dtdt 

^ti . (5a) 

Th 6 in'bcgral on "bho right can he evaluated since the integrand and limits 
are known. For brevity, it will be designated Q 12 . 

Z2 = zi + wi (ta - ti) + — ^12 (51>) 

In equation (5t>) “the unknowns are and W i ; the_^conjtgnt Wi — elim- 

inated by using data from a third station: 

Z 3 = zx + wi (ta - ti) + Qx3 (5c) 

solving equations (5b) and ( 5 c) simultaneously yields the following estpres 
Sion for 




Zg - Zi - (Z3 - Zi) . ti 


Qi2 - Q13 


tg - tx 
ta - ti 


(6) 


In order to determine the center-of-pressure position, it was also 
necessary to measure the pitching-moment-ctirve slope, Cma,, for use in the 
equation relating lift- and moment-curve slopes: 


1 M ct - C™ 

I '-Lq, - '-mo, 


( 7 ) 
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The value was obtained from .the frequency of oscillation as 

determined hy the least-squares fit to the angle -of -attack history. The 
equation for the frequency of oscillation with one degree of freedom was 
used: 


^ _ 1 _ 1 y-Cma,<loSZ 

^ 2n^ ~ 2«V I 

or 


^ ~ * Sqo^ 


(8) 


The method of calculation of lift requires data from three stations 
in the wind tunnel to determine one lift answer. For each of the tests 
from which lift-curve slope was obtained in this investigation, four shad- 
owgraphs were made. This permitted making four combinations of three 
stations from each set. Occasionally, one combination of three stations 
failed to define the lift-curve slope because the position of the stations 
along the flight path was such as to record three values of z which fell 
nearly on a straight line. Such badly defined values are not included. 
Typical data from which lift and center-of -pressure location were deter- 
mined are shown in figure 1)-. • 

Thus far the motion has been treated as an oscillation in the x,z 
plane. Actually, the models were free to move in the x,y plane as well. 
If the aerodynamic properties of the model were perfectly linear, there 
would be no influence on the data if this fact were ignored. Since it 
was expected, in this case, that the lift would be influenced by sideslip 
single, the value of Ci^ was related to the mean absolute angle of yaw, 
the vector sum of a and 3, diiring the interval in question, rather than 
to a alone. Thus, 


where is called the lift coefficient corresponding to the mean 

absolute angle of yaw, is the value of lift-curve slope obtained as 

described above, and ap is the mean absolute angle of yaw for the test 
interval . 


ACCURACY 


A realistic estimate of the accuracy of the above techniques must 
of necessity be based on consistency of results, rather than on a detailed 
study of all factors contributing inaccuracies. The lift results are 
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consistent within a scatter in Clr of ±0.008, based on frontal area, 
for the present test. This amounts to less than ±4 percent of the meas- 
\ired lift at 5° angle of attack. The drag is believed accurate within 
2 percent. The total scatter in center-of -pressure position of six sep- 
arate tests was ±1.2 percent of the body length. The static stability 
margin of the models was made small, from 0.3 to 2.5 percent of 1; 
therefore the determination of center-of-pressure position was not sen- 
sitive even to large errors in Chiq^. 

RESULTS 


The testing technique is such that varying amoiints of disturbance 
are imparted to the models as they are launched. These differences, plus 
slight damage suffered by a few models, made it necessary to consider 
each test individually. Only the aerodynamic parameters best defined by 
each test have been included here as explained below. 


Lift 


The variation of lift with angle of attack, Clr versus aj^, as deduced 
from four tests (l, 2, 3, and 5 )>_ls presented, in. figure-5. — Out— of-a 
total of " l6 “lift -curve -slope determinations, 2 were deleted, 1 each from 
tests 2 and 5, because the data were such as to fail in defining the 
flight -path curvature accurately. 

Test 1 had a peak amplitude of oscillation of about 1° and deviated 
about 0.07 inch from straight-line flight. Tests 2, 3> and 5 had larger 
amplitudes. The model in test 2 was bent about 1.5° at the base of the 
spike during the launch (see fig. 6(b)). The effect of this damage was 
estimated by calculating the theoretical lift on the conical spike nose 
(ref. 3) at 1.5° angle of attack and doubling this to allow for lift 
carry-over on the. spike cylinder. This lift was then added in coefficient 
form to the values of Clr determined from this test and amounted to a 
correction of about 7 percent of the final answer. Test 3 had the largest 
amplitude recorded in the investigation - about 10°. The nose tip (about 
1 percent of the model length) was hooked over by lateral acceleration in 
the gun. The hook was bent in a plane nearly normal to the x,z plane, 
so it is believed that this damage had negligible effect on the lift. 


Center of Pressure 

Tests of five models (l, 5, 7, Q, and 9) showed the center of pres- 
s\ire to be between 72.1 and 74.6 percent of the length aft of the nose. 
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For all of these tests the oscillatory motion was fairly well defined 
and the models were not visibly damaged (see fig. 6). 


Drag 


Five tests were obtained with accurate enough time-distance measure- 
ments for determining drag (l, 4, 5, 6 , and 7). These results are plot- 
ted in fig\are 7 as a function of the average position of boundary- layer 
transition along the model as observed in the shadowgraphs. The drag 
due to lift of these models was estimated and subtracted as outlined in 
reference 1. 

The principal evidences of txirbulent boundary- layer flow are eddies 
visible in shadowgraphs and small waves sent out by these eddies. A region 
of tiarbulent flow may be seen near the base of the model in figure 6(j). 
This patch of turbulent flow is being swept downstream at about half the 
free-stream velocity as evidenced by the angle of the weak shock wave 
emanating from it. Scrutiny of all the pictures from a test permitted 
the estimation of transition position. In the case of test 6 no laminar 
flow was seen aft of the spike cylinder. It is believed that transition 
did, in fact, occur at the nose due to some slight damage not visible in 
the shadowgraphs. In view of the difficulty of seeing boundary- layer 
turbulence near the nose, the point for this test. is plotted as if tran- 
sition occurred at the spike base. 

The drag' tests of two unmodified cones are also reported in this 
figure. The lower Reynolds number of these two tests resulted from using 
shorter models and testing at slightly lower tunnel pressure for testing 
convenience. 


DISCUSSION 


The initial lift-curve slope predicted using reference 3 is indicated 
in figure 5 and is seen to iinderestimate the lift appreciably at angles 
of attack greater than 3°. A somewhat better model of the flow at angle 
of attack is given by impact theory (ref. 4) which predicts increasing 
lift-curve slope with increasing angle of attack. This predicted varia- 
tion is plotted in figiire 5; there is still a large gap between ej 5 )eriment 
and theory. 

The nonlinearities indicated in impact theory may be evaluated dif- 
ferently using the crossflow concept proposed in reference 5 and further 
developed in reference 6. If the crossflow contribution is added to the 
continuation of the initial lift-curve-slope line of figure 5, the sum 
agrees quite well with the experiment. This crossflow- lift increment was 
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calculated on the basis of TO percent of the supercritical drag coeffi- 
cient of a circular cylinder being developed on an expanding body at 
angles of attack above 2.5°, This value, 70 percent, was chosen on the 
basis of results given in reference 6 for the nose of an ogive cylinder. 

The agreement, perhaps better than might reasonably be expected, still 
indicates the possibility of accurate lift calculation by presently avail- 
able methods . 

The effect of the spike on the lift appears to be small. Its frontal 
area, less than k percent of the model base area, would indicate that even 
if the lift carry-over from the nose were equal to the nose lift, the 
total lift would be increased by only k percent. The spike apparently did 
not affect the lift on the main body to any great degree, as was indicated 
by the success in calculating the total lift. 

The center-of -pressure position was estimated theoretically, ignoring 
the lift developed by the spike cylinder and assuming that the pressures 
on the aftercone were not affected by the presence of the cylinder. The 
resulting center of pressure, essentially the center of plan-form area of 
the conical portions, is predicted to be 7k. 2 percent of I aft of the 
nose. The median experimental value of 73.3 percent agrees well with this 
prediction. If the lift carry-over from the nose cone onto the cylinder 
is assijmed equal to the nose-cone lift and acts at one third the cylinder 
length aft of the nose cone, the predicted location is 73* k percent I 
aft of the nose. It is expected that at large angles of atta ck..— the center 
of pressure -moves— forward'”slightly since the spike-cylinder lift variation 
is expected to be more nonlinear than that for the conical portions. 

The estimates of the zero-lift drag of the projectile were based on 
references k, 7, and 8 for wave drag and reference 9 for skin friction. 

The base pressure was assumed to be 10 percent of free-stream static 
pressure, based on extrapolation of the data of reference 10. The wave 
drag estimated by use of reference k was about 20 percent below that pre- 
dicted by use of reference 7 for a 5 ° cone. The value from reference 7 
only is presented in f igvire 7 • 

The theoretical drag build-up for both config\irations is presented 
in figure 7. The skin friction was calculated using local conditions and 
the effect of heat transfer was included. The data of reference 9, and 
subsequent data, as yet unpublished, obtained by the same technique, were 
used in these calculations. 

The effect of heat transfer was to increase the skin friction by 
some 35 percent over that which would have existed without heat flow, 
according to reference 9* In^ order to see if the wall to free-stream 
ten^jerature ratio was the same as for the tests of reference 9» estimates 
were made of the maximum temperature rise the model could have experienced. 
Before firing, the model ten^serature was around 5^° R "the free-stream 
static temperatiire was 300° R. At M = 8.6 the recovery ten 5 >erature was 
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around 1^300° R, and the flight duration was 0.008 second. At the outset, 
it is useful to note that even if the entire model surface were to be 
raised to the melting point, iteO® R, the temperature difference causing 
heat transfer would be reduced by only 23 percent. That this extreme 
temperature rise did not, in fact, occur is evidenced by the perfect con- 
dition of models observed at the last shadowgraph station. 

Calculations of the maximum possible temperature rise of the model 
based on the appendix of reference 11 indicated a temperature rise on the 
spike cylinder, just behind the nose cone, of 30° R and 1000° R at the 
apex. Thus the temperature rise on 97 percent of the wetted area was 
less than 30° R and the maximiam temperature rise at the nose, as indicated 
by its survival, was less than 880° R. Thus it is seen that the wall to 
stream temperature relation of the skin-friction tests also applies to 
the present test. 

Two calculations of the local Mach number and dynamic pressure, on 
which the skin-friction calculations were based, were made. The first 
was based on reference 7 and ignored the effect of the spike on the flow 
over the main body. The second was based on reference 7 for the nose cone, 
and shock-expansion theory, reference 8, for the remainder of the body. 

The pressxire distributions and skin-friction results from these two cal- 
culations were nearly identical, indicating no important effect of the 
cylinder on the pressure distribution on the remainder of the body at 
small incidence. 


As is indicated in figure 7, boimdary- layer transition generally 
occurred well downstream of the spike base. In many cases the flow was 
laminar over the entire body length, but it never remained so throughout 
an entire flight. Thus it appears that at the present test conditions 
of high Mach number and large heat-transfer rate into the model surface, 
a laminar boundary layer is extremely stable. 


CONCLUDING REMARKS 


Presently available theories may be used for engineering estimates 
of lift, center -of -pressinre location, and wave drag of the projectile and 
the drag of the 5° half -apex-angle cone from which the projectile shape 
was derived. The lift for the projectile at angles of attack around 3° 
and above is noticeably greater than predicted by linear theory. 

A good estimate of friction drag was made on the basis of presently 
available data which include the effect of heat transfer. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 23, 195^ 
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Figure 1. - Models and sabot. 
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Figure 2.- Expanded assembly of model and sabot. 



A-18381 

Figure 3»- Photomicrographs of model surface and profile, XUOO. 


CONFIDENTIAL 










Present test 



ONFIDEKTIAL 


Figure 5 . - Variation of mean lift coefficient with mean angle of attack. 
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(a) Test I , station 2 



(c) Test 3, station 4 



(d) Test 4, station I 



(e) Test 5, station 2 


Figure 6.- Shadowgraphs 



(f) Test 6 , station 3 



(g) Test 7 , station 4 



( i) Test 9 , station 3 



(j) Test 10, station I 



(k) Test II, station I 
of models in flight. 
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Figure 7.- Variation of drag of projectile and a 5® cone 
with position of transition location. M 8.6. 

COHFIDEHTIAL 


NACA-Langley - 11-19-54 - 425 




